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We present a Random-Phase- Approximation formalism for the collective spectrum of two hyperfine 
species of dilute 40 K atoms, magnetically trapped at zero temperature and subjected to a repulsive s- 
wave interaction between atoms with different spin projections. We examine the density-like and the 
spin-like oscillation spectra, as well as the transition density profiles created by external multipolar 
fields. The zero sound spectrum is always fragmented and the density and spin channels become 
clearly distinguishable if the trapping potentials acting on the species are identical. Although this 
distinction is lost when these confining fields are different, at selected excitation frequencies the 
transition densities may display the signature of the channel. 
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I. INTRODUCTION 

The realization of Bose-Einstein condensation of alkali 
atoms in magnetic traps triggered substantial experimen- 
tal efforts, aimed at cooling fermion isotopes below their 
Fermi temperatures. In particular, lithium and potas- 
sium have been trapped and cooled Jj]-^, and a recent 
experiment displays unambiguous evidence of quantum 
degeneracy effects at temperatures around half the Fermi 
temperature of 40 K atoms Q in a harmonic trap. On the 
other hand, various theoretical descriptions of thermo- 
dynamical properties of confined, cool free fermions have 
been presented, either in the semiclassical limit Jp) or 
with explicit consideration of quantum shell effects pf . In 
Ref. a Hartree-Fock (HF) calculation of the mean field 
spectrum of two hyperfine fermionic species subjected to 
s-wave attraction has been forwarded. Moreover, in view 
of the presence of both fermion and boson isotopes in 
natural alkali samples, the consequences of their coex- 
istence and mutual coupling at zero temperature in the 
magnetic trap is a topic of current interest |§|-|l3| , as well 
as the possible occurrence of BCS-like superfluid states 
driven by attractive interactions [p4|-p~9| . 

An important step towards a full understanding of the 
behavior of coexisting hyperfine species is the knowledge 
of their collective excitation spectrum. In this context, 
the collisionless modes of an extended system with vari- 
ous hyperfine levels have been examined in the frame of 
Landau's theory of Fermi liquids p0| and the zero-sound 
collective spectra of two species of confined fermions 
have been computed resorting to sum rules |21|. As- 
suming local equilibrium of a Fermi gas, described by 
a Thomas-Fermi (TF) approximation, the linearized hy- 
drodynamic equations can be analytically solved both in 
the degenerate and in the classical limits |2^| . A related 
hydrodynamic-like approach based on the equations of 
motion for the first and second moments of the fermionic 
Wigner distribution allowed to compute the oscillation 
modes of an isotope with one E3] and two spin compo- 



nents |2J. 

Since the two hyperfine species of trapped fermionic 
atoms are very dilute, it is not clear that when the system 
is excited by a low frequency external field, the s-wave 
interaction -which is supposed to play a relevant role in 
thermalization during the evaporative cooling process- 
may permit the trapped gas to equilibrate locally and 
develop hydrodynamic oscillation modes. In particular, 
it has been shown that for 40 K, at least 10 8 atoms should 
be cooled in each hyperfine state in order to reach the 
hydrodynamic regime in the degenerate quantum phase. 
It is then important to focus upon the study of colli- 
sionless excitation spectra of these systems, seeking to 
improve the understanding of their low temperature be- 
havior, as well as the evolution between the collisionless 
and the hydrodynamic regimes, as increasing tempera- 
ture suppresses Pauli blocking effects and enhances the 
collision rate. For this sake, in this work we derive a 
Random-Phase- Approximation (RPA) description of the 
collective modes of two species of fermions in a harmonic 
well with mutual s-wave coupling at zero temperature, 
and apply the formalism to the computation of density- 
like and spin-like fluctuations. The paper is organized 
as follows: the specific RPA frame and the extraction of 
the elementary excitation spectrum of quasiparticles in a 
mean field approach are discussed in Sec. [n . In Section 
III we propose a simple model to describe the zero sound 
modes that permits to analyze various major features of 
the expected RPA spectra. Typical calculations of collec- 
tive spectra for the lowest multipolarities are presented 
and discussed in Sec. IV, while Sec. 
main conclusions. 



summarizes our 



II. THE RANDOM-PHASE-APPROXIMATION 
FOR A TRAPPED FERMION SYSTEM 

We assume that the trapped atom system consists 
of nonintcracting quasiparticles (qp's) in a mean field. 



I 



Throughout this paper this is referred to as free sys- 
tem, which can be excited by an external field so that 
particle-hole (ph) pairs involving, in principle, both hy- 
pcrfine species (oi,^), are created with energy SI. The 
spectral properties of this nonhomogeneous free system 
are contained in the free ph propagator Gq a (SI), where 
the labels a and a' stand for either o\ or a 2 . 

We also suppose that a ph effective interaction 
Vph a a acting between the qp's gives rise to the 
dressed propagator for ph pairs according to the RPA 
integral equation |p5[ 



G° 



(O) = Gq «' (Si) 



G T 



\si) (2.i) 



It is worthwhile noticing that the system under con- 
sideration is both nonhomogeneous 26-28 and polarized 
p9| , ^p| , so that the present development merges the corre- 
sponding formalisms as shown below. Hereafter, we con- 
sider longitudinal excitations involving propagation of ph 
pairs of the same spin kind, created by spin-symmetric 
(s) and spin-antisymmetric (a) multipolar operators of 
the form 



2Vi 



N 2 
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with 
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(2.2) 
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being N\ (resp. iVjj) the number of trapped atoms of 
species o\ (resp. 0-2)7 an d L the multipolarity of the 
perturbation. Notice that a particle-particle interaction 
V(r — r'), with r, r' respectively denoting particles with 
spin projections a and a', gives rise to a ph interaction 
that scatters a (aa^ 1 ) pair at position r into a (a'a'~ ) 
pair at r', and that only collisions among different species 
are allowed . The free ph propagators involved in longi- 
tudinal density fluctuations are diagonal in spin space 
and thus Eq. (2.1) splits into two equivalent systems of 
two equations each, intrinsic to polarized systems |^9|,^0| , 
which in coordinate representation read 

G^ 1 (r,r') = Gr"(r,r') 

+ J d 3 ri d 3 r 2 G^'"(r,r 1 )V/ I Z'(r 1 ,r 2 )G CTV "(r 2 ,r') 
G CT ' CT " 1 (r,r / ) = 

d 3 ri d 3 r 2 G£ V " (r, n) vX'in, ra) (r 2 , r') 



(2.4) 

In spatially inhomogeneous systems, it is convenient to 
expand both free and dressed propagators in multipolar 
decompositions Bqjgjl 



G ctV " (r, r', SI) = £ Gl a " (r, r', SI) P L (f ■ f') (2.5) 

L 

with the Legendre polynomials Pl (x) . The free ph prop- 
agator reads 

Gr"(r,r',0) = 

£ €(r) [€>(rT mr')]* XZAV) (2-6) 

vv' 

with (f>° (r) a single-particle (sp) wave function for energy 
eigenvalue and xZu> ^ ne generalized susceptibilities in 
terms of the Fermi-Dirac occupation numbers n(s) 



n{el)-n{el,) 
n-(e£-eZ)+*rj 



(2.7) 



Here the label v stands for the spherical quantum num- 
bers (nlm). The expressions for the multipole compo- 
nents Go J* (r, r', SI) are given in the Appendix. 

Since (2.4) is a matrix equation calling for discretiza- 
tion in radial coordinates, it is convenient to map it onto 
a vector system for the transition densities defined as 



ill' (r, SI) = J dr' G°'°~ X (r, r', SI) r' L Y LM (f') 
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where 
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The multipolar susceptibility can then be computed as 
' ~\SI)= f dvr L YU{f)Spi^\T, Si) 
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For dilute trapped systems at low temperature we 
can reasonably represent the interaction potential by a 
contact interaction of the form gb(v — r'), being g = 
AttTl 2 a/m with m the mass and a the s-wave scattering 
l eng th of the interacting atoms. Thus, we obtain from 

8pl"~ 1 (r,Q) = 8fi5f 1 (r,Q) 
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(r,r\Sl)5pl°- (r\Sl) (2.11) 



Moreover, in view of (2.2) we shall consider the sym- 
metric and antisymmetric density fluctuations for each 
atom species, 
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This representation permits to distinguish the here- 
after called density-like (symmetric) and spin-like (an- 
tisymmetric) fluctuations |pl| as they usually appear in 
Fermi liquids. In fact, a multipolar operator 0\ will likely 
generate a total density-like fluctuation proportional to 
5p s a + 8p s a , while out-of-phase perturbations produced 
by 0\ will induce a spin-like fluctuation proportional to 
&p% + 8p a a , . In addition, this distinction enables one to 
analyze the influence of a given fluctuation on its own 
propagation, as well as on the oscillations in the other 
spin species. 

The numerical pr oced ure consists of solving the dis- 
cretized equations ( 2.1l|) by matrix inversion, comput- 
ing the susceptibilities (P-10| ) and constructing the total 
dynamic structure factors S*( s ^ a ) = — Imx^ s,a ^/7r in both 
spin channels M,B(|, where 
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(2.13) 



The collective spectrum of density-like and spin-like 
modes for a given multipolarity L is indicated by the 
poles of the real part of these responses, or correspond- 
ing peaks in the dynamical structure factors. 

It is important to remark that this is a very general 
RPA description of collective excitations, valid for any 
system identified by its elementary excitations with sp 
spectrum e„ and states <j> v , and by an effective ph inter- 
action Vph- I n most applications to quantum liquids (see 
for example Refs. p7|-p0[| and therein), one starts from 
a HF sp eigenspectrum and chooses the ph coupling as 
the double functional derivative of the total energy with 
respect to the sp density. The HF spectrum of two hy- 
pcrfine species of trapped fermions has been previously 
investigated in Ref. M for the case of an attractive cou- 
pling between the species, and in the present work we 
adopt the same philosophy for a repulsive interaction of 
strength g. The HF spectrum arises from the solution of 
the coupled nonlinear system in spherical coordinates 



2m dr 2 
(r) 



h 2 1(1 + 1) 
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for species a ^ er', with partial densities 

\< t (r)\ 2 



nl 
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(2.14) 



(2.15) 



and trapping potentials fnui 2 r 2 /2. 

For vanishing temperature, the Fermi-Dirac occupa- 
tion numbers are step functions limiting the summation 
to states below the respective Fermi sea e F that fulfills 
the number equation N a — ^ v <d(sp — e a v ). In the forth- 
coming calculations, we shall consider both equal as well 
as different trapping frequencies uj a for species 1 and 2, 



the latter case devised to take into account the corre- 
sponding magnetic projections of the trapped atoms, i.e., 
K/^ 2 ) 2 =ai/a 2 . H 



III. A SIMPLIFIED MODEL FOR ZERO SOUND 
MODES 

In order to get a basic understanding of the differ- 
ent excitations of the two-component gas, we propose a 
very simple model for the noninteracting propagator that 
makes room to illustrative analytical results. Let us as- 
sume that each of the coexisting species possesses only 
one elementary excitation with energy u a (resp. clv). 
The corresponding free propagator reads 



Gr(r,r') = F CT (r) [F<V)]* xZ 



with 



1 
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(3.1) 



(3.2) 



being ^(r) the wave function of the exc ited ph pair at 
the given position. Rep lacement of ( |3.l| ) into the RPA 
system of equations (2.4) brings into evidence that these 



propagators are of the form 

g™-\t,t') = f°{t)xZ [r'M 

G-'°- 1 (v,ri) = gF°'(r) X a Q ' [^''"V) 
where in turn 



(3.3) 
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(3.4) 



represent dressed ph wave functions, which depend on 
the ph interaction strength 



,/ = // / dx[F-(x)]* F° (x) 
(aa-^Wa'- 1 ) 



(3.5) 



The system (3.4) possesses a simple algebraic solution, 
from which we can write the transition densities as 



8p°°- (r) 
5/" -1 (r) 



Xg [F a (r)] * 8° 
D 

- Xlxl [F°(r)}* 0° 
9 D 



being 



£>M = l-|5| 2 XoHXo (w) 



(3.6) 



(3.7) 



with Q a = J d,Vi 0\ F a (r i ) the matrix element of the tran- 
sition operator. Consequently, the dynamical susceptibil- 
ities are 
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The appearance of the common denominator D{u>) in 
Eqs. ( |3.§| ) reflects the fact that the RPA describes collec- 
tive fl uctu ations of the system as a whole. The zeroes of 
D in (3/7) can be easily found to be 



± 



16\g\ 2 u a 



1/2 



(3.9) 

The evolution of the collective modes with interaction 
strength is encompassed in Eq. ( |3.9| ); assuming for in- 
stance uj 2 > u)p, , to lowest order in the expansion pa- 
rameter one finds, 



0J„ 



2 2 



A\g\ 2 uj a uj a 



A\g\ 2 uj a u a 
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(for real g). It is i mportant to remark that although both 
frequencies (3.13) correspond to poles of each term in 
( |3.14 ), explicit computation shows that w CT + \g\ has zero 
amplitude in the antisymmetric response while uj„ — \g\ 
has zero amplitude in the symmetric one. The higher and 
lower frequencies in ( 3.13 ) can be then attributed to the 
symmetric and antisymmetric fluctuations respectively. 
Their intensities, at the corresponding energies, can be 
derived as 



4 ^ 
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The latter line holds in view of ( 3.12 ) and coincides with 
the unperturbed strength of the elementary excitation. 

In summary, this simple model gives rise to s pect rum 
fragmentation into two lines, according to Eq. ( |3 . 9[ ) ; in 
the specific case of equal trapping frequencies and in the 
weak coupling limit, we see that due to amplitude can- 
cellation, one can clearly identify different excitations en- 
ergies in the s and a channels. 



showing that within the validity of the approximation 
involved, the collective modes progress towards opposite 
directions with increasing coupling. For large interaction 
strength and/or close intrinsic frequencies of the sp ecie s, 
one can derive a complementary limit (cf. Eq. ( [3.9] )) 
which gives 



± 2\g\y/u} a uj 



32|.g 



n\2, 



(3.11) 



The particular situation where the two interacting 
species are identified by the same well parameters is es- 
pecially enlightening. On the one hand, one can see from 
(3.11) that if uj a = Lo a i, the dispersion relation of the 
modes is simply given by 



1 ± 



(3.12) 



which for sufficiently small relative strength \g\/uj a illus- 
trates the fragmentation of the elementary excitation into 
two collective states with energies 



IV. CALCULATIONS AND RESULTS 

We have solved the HF problem for variable numbers 
of 40 K atoms N\, N2 corresponding to different spin pro- 
jections in an isotropic harmonic trap, with a mutual 
s-wave scattering length a — 8.31 nm |3[. As in Ref. [0, 
we start an iterative procedure from oscillator wave func- 
tions «°f (r), and convergence is rapidly achieved. The 
sclfconsistent states are labelled by the same quantum 
numbers nl and the wave functions differ only slightly 
from the original ones. The low energy states are the 
most sensitive to the size of the interaction strength, 
reaching deviations with respect to bare oscillator en- 
ergies as large as 15%. The combinations of ph states 
entering the free ph propagator (see Eq. ( |A3|) ) are se- 
lected by angular momentum conservation, and we find 
that the elementary excitation energies in the denomi- 
nators of the generalized susceptibilities (see Eq. (|2.7|)) 
are weakly spread around the noninteracting oscillator 
values. 



A. Equal trapping potentials and populations 



UJ = LO a ± 



(3.13) 



The corresponding form of the total spin symmetric and 
antisymmetric responses is (cf. Eq. (2.13)) 



/(»!<») _ 



4 ^ 2(1 ±3x0) 



(3.14) 



As a first step we examine a trapped two-component 
Fermi gas with equal trapping potentials (ETP), i.e. 
u)\ = u>2, and the same number of atoms in each hy- 
pcrfine level. Under these conditions the RPA equations 



Q2.4|) can be decoupled for Spa ^ and Sp 1 ^' , giving rise to 
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(s,o) 



Spi 



and Gq l 



Ga a 
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for either spin projection a. Equations (4.1) are iden- 
tical to the RPA equations for a single component gas 
interacting through a repulsive (resp. attractive) contact 
interaction, giving rise to the symmetric (resp. antisym- 
metric) fluctuation. 

In Fig. |l| we show the strengths S^'^ of the monopolar 
(L = 0) fluctuations for N12 = 10 4 . It should be kept 
in mind that the spread of the peaks in S^ s ' a ^ (Q) is an 
artifact of the calculation, due to the in trod uction of the 
small numerical parameter r\ (cf. Eq. (2.7)). The large 



scale behavior of S(Q) shows multiple excitations at en- 



ergies close to ft r , 



e n2 i w 2nui . These peaks 



decrease their amplitude (notice the logarithmic scale) as 
we increase the transferred energy. Although the s and a 
channels are different, we are not able to visualize them 
in the current scale. A narrow region around the oscil- 
lator excitation energy is displayed in the lower plot of 
Fig. [I]. In fact, a careful analysis of the peaks indicates 
that they group into two main sets. Only the strongly 
fragmented one around £1 = 2lj\ can be viewed in this 
figure; the intensities in the second group lying at higher 
energies are too small in the current scale. Moreover, the 
spectra in both channels arc clearly different: while the 
s channel is shifted upwards in energy, the a one lies at 
lower energies, in qualitative agreement with the simple 
model estimate of the preceding section. It is important 
to mention that the qp's energies coming out from the 
HF calculation are such that the ph excitation frequen- 
cies are lowered with respect to the oscillator ones by the 
mean field interaction, while the ph coupling introduced 
in the RPA formalism shifts the collective spectrum to 
higher energies. In addition, we have numerically veri- 
fied that the position of the main peak can be accounted 
for by the analytical sum rule formula derived by Vichi 
and Stringari pi| . 

The behavior of the dipolar fluctuation (L = 1) is 
slightly different. Although both the s and a spectra 
are weakly fragmented, the amplitude of the oscillator 
mode at fi = u)\ is large. It is well known j}3) that for 
equal trapping potentials and populations there is a dipo- 
lar excitation in the symmetric channel, associated to the 
center-of-mass oscillation of the gas and occurring at the 
oscillator frequency. However, in the spin-like channel 
this excitation appears at a slightly lower energy. These 
facts arc verified in the present RPA calculation as shown 
in Fig. |2j where we see the dipolar structure factor for 
N\ = N2 = 10 4 atoms. We can also observe that in 
addition to these modes, fragmented poles of collective 
nature show up as indicated by the simple model of the 
preceding section. 

Let us now examine the spatial profiles of the den- 
sity fluctuations. In principle, one would expect that if 
the s yste m is excited with a given operator s or a (cf. 
Eq. (2.2)) in the vicinity of a peak in one of these chan- 



nels, a density fluctuation will develop that reflects both 
the character of the external field and the nature of the 
intrinsic excitation of the free system. As an illustration, 
in Fig. H we show the real parts of the monopolar density 
oscillations for the same conditions as in Fig. [l] at the 
frequencies given by the poles of the response. We ob- 
serve that when the symmetric fluctuation is important 
within the bulk of the trapped system, the antisymmetric 
counterpart is completely negligible and viceversa. This 
is in agreement with the interpretation of the simplified 
model of Sec. Ill in which density and spin-like fluctua- 



tions were clearly distinguished. However, at the energy 
of the pole in the symmetric channel (f2 s « 2.03), surface 
oscillations develop both in 8p s and 5p a ; this behavior 
can be attributed to the proximity of a weak peak in the 
antisymmetric structure factor (not visible in the scale of 
Fig. 1). 



B. Unequal potential wells and populations 

In the recent experimental achievement of DeMarco 
and Jin 0], 40 K atoms in two magnetic sublevels \F = 
9/2, m F = 7/2) (type 1 atoms), and \F = 9/2, m F = 
9/2) (type 2 atoms), were confined and cooled. Although 
the ratio oji/ui = (02 /ai) 1 / 2 ~ 1.13 is very close to 
unity, we take it into account as an explicit feature of 
the real nonsymmctric configuration. We also consider 
unequal populations ranging from A = N2/N1 0.3 to 
A w 3, which can be built during the evaporative cooling 
process. 

As in the ETP case we analyze the dynamic struc- 
ture factors and the fluctuations created by multipolar 
external fields. In Fig. || we show the symmetric dy- 
namic structure factor for several values of the popula- 
tion ratio A, N\ = 10 4 and a monopolar excitation. The 
structure of spin-like fluctuations cannot be distinguished 
from the density- like ones within the scale of this plot. 
The essential features of the response can be summarized 
as follows: i) in either channel, the strength of each ph 
transition is redistributed so that the structure factor is 
fragmented around the noninteracting peaks with some 
intensity appearing at higher energies; ii) for low or high 
values of A the structure of the most populated species 
essentially reproduces the pattern of the noninteracting 
quasiparticles, while the spectrum of the other species is 
highly fragmented with a largely suppressed amplitude, 
iii) given either species, as the number of atoms in the 
other spin projection increases, the complementary exci- 
tation appears displaying considerable fragmentation. 

We have also made calculations for multipolar excita- 
tions with L = 1, 2 ; in Figs. || and || we show the corre- 
sponding results for the structure factor. In general, the 
description of the monopolar excitation applies as well 
to higher multipolarities, however we can mention some 
differences: the dipolar peaks are narrower and fragmen- 
tation seems to be stronger. In particular, in the equal 
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population case (A = 1) there exist well resolved frag- 
mented peaks as intense as the original HF ones. 

The real parts of the transition densities are displayed 
in Figs. to H for L = to 2, respectively. The s and a 
density profiles for each component have been scaled to 
a common value at a given frequency. This is a conve- 
nient criterion, since if one depicts the four profiles in the 
same scale, near the intrinsic ph frequencies of each com- 
ponent the complementary density fluctuation appears 
largely depleted. Before analyzing the fluctuation pro- 
files, we want to bring some attention to the shape of the 
HF equilibrium density profiles. In Table | we quote, as 
a function of A, some shape parameters related to the 
spherical probability density distribution P a (r) defined 
as: 



r 2 p a {r) 



(4.2) 



Particularly, we list the maximum probability P max 
which is attained at i? ma x, the location of the density 
edge i? c dgc and the full width at half maximum (FWHM) 
of the probability distribution. These parameters will 
help us to understand the main aspects of the excitation 
profiles. 

In Fig. [7] we show some typical monopolar density fluc- 
tuations for each spin component, at the lowest-lying col- 
lective peaks and several concentrations. We observe that 
as we increase A (from bottom to top in Fig. Q), both 

fluctuations 5p 2 s ' a ^ for a given frequency extend beyond 
the P2 edge. As an illustration, let us consider in detail 
the case A = 0.3. For the excitation with SI w 2wi, 
8p2^ is bounded to roughly the size of the density of er 2 

atoms (ar < 7, a 2 = h/muji ) while 6p^' a ' extends far 
beyond this cutoff, within the p\ range. For small A the 
corrections Sp 12 are unimportant, in general, due to 
the relatively weak coupling of the <j\ species to the few 
02 atoms. However, in the same limit Sp 22 acquires 
important corrections 5p 21 , essentially driven by the 
large number of <j\ atoms. We also see two different be- 
haviors related to the SI dependence: for frequencies close 
to 2u>i, 5p2 and 5p 2 present opposite signs, revealing 
that the induced contribution 5p 2 1 is larger than the 
intrinsic one Sp 2 2 . In fact, if the system is excited with 
frequencies close to 2 lo\ , large Spi fluctuations should be 
expected, which in turn introduce, through the V p h inter- 
action, sizeable Sp 21 contributions to Sp 2 s ' a ^- In turn, 
exciting at SI w 2ui2 creates Sp 22 \ thus the correction 
to Sp 2 1 is second order in the interaction and the s and 
a fluctuations are the same. On the other hand, as p\ ex- 
tends beyond pi , a stimulus acting at SI = 2 0J2 induces 
a cross-fluctuation 5p 12 inside the bulk of the type 1 
system, i.e., at smaller radii (cf. Fig. [?]). 

As one increases the number of particles in species 2, 
their density profiles extend farther (see Table |) and 
there are few noticeable changes in the spatial localiza- 
tion of the excitation. However, in general, close to 2u>2 



we find a strong induced fluctuation 8p 12 in the spatial 
region where P2 (r) is larger than Pi (r) and correspond- 
ingly for Sp 21 1 near 2u>i. If one keeps increasing A, the 
perturbative interpretation is no longer valid for every 
multipole and species. For A = 3, the fluctuations oc- 
cupy a broader region, roughly from ar = 4 to ar = 10, 
this can be attributed to a wider P2 density. 

In addition to the vanishing at origin of the L 7^ fluc- 
tuation, the main difference between the spatial profiles 
of distinct multipolarities, displayed in Figs. [?] to |], lies 
in the enhanced symmetry of the channels for L = 2. In 
this case, we observe that the cross-fluctuation of a given 
component is non-negligible only near the free ph transi- 
tion of its counterpart, suggesting a weaker p h effe ctive 
interaction for higher multipolarities (cf. Eq. (2.11)). 

Another feature to remark is the difference with the 
ETP case; in that situation we have the same spatial 
fluctuations for either spin projection and well defined 
channels s and a; however, in this more general prob- 
lem, we observe different spatial profiles for each compo- 
nent and mixed behaviors within the channels, i.e., given 
a pole in the symmetric susceptibility, fluctuations may 
exhibit similar amplitudes in both channel. 



V. DISCUSSION AND SUMMARY 

In this work we have developed a RPA formalism for a 
two-species, trapped Fermi gas at vanishing temperature, 
which may provide guidelines to current experimental re- 
search. We have shown that the interspecies interaction 
gives rise to a fragmented zero sound spectrum. We ana- 
lyzed both the equal trapping potential and populations, 
as well as a general case with unequal potentials and pop- 
ulations intended to mimic an experimental situation. 
The main differences arise in the spectra; while in the 
ETP case the poles associated to density-like and spin- 
like fluctuations were clearly distinct, in the general case 
these fluctuations cannot be disentangled. The density- 
like and spin-like responses of the system share the same 
energy spectrum, however with unequal amplitudes. We 
have presented a simplified picture for the ph propaga- 
tor and elementary spectrum which gives support to the 
numerical results. This mo del al so allows us to interpret 
the frequency shifts in Eq. (3.13) in terms of the interac- 
tion among pair excitations, measured by the product g 
of the bare coupling strength time s th e overlap between 
the pair wave- functions (cf. Eq. (3.5)). For the excita- 
tion with L = 1 we have found in the density channel 
and for the ETP case a pole at the bare oscillator fre- 
quency corresponding to rigid oscillation of the system; 
however this is no longer true for a more realistic con- 
figuration with unequal trapping frequencies . Al though 
we checked that the eigensolutions of Eq. (2.14) differ 
only in a few percent from the bare oscillator basis func- 
tions, we cannot prevent propagation of small errors in 
the response calculation, which ought to be safely com- 



G 



puted out of the true HF wave functions to the price of 
a increase in computing time. Finally, we should men- 
tion that this two-component RPA formalism may be 
straightforwardly generalized to nonzero temperatures by 
taking into account the finite temperature Fermi-Dirac 
occupation numbers. 
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APPENDIX: 

To evaluate the free ph propagator we use the 3D har- 
monic oscillator basis for each spin projection: 

#UmW = Kl exp(-a 2 r 2 /2) r l L^\a\ r 2 ) Y lm (r) 
= Rni(r)Y lm (f) (Al) 

being a 2 = h/muja-, Yi m spherical harmonic functions, 
_Lri +1//2 the generalized Laguerre functions and A a nl a nor- 
malization constant. The multipolar component Gq£ 
reads 



Gil (r,r',Q) = 



1 



(47T) 



(2l+l)(2l' + l) 



n l.n' V 



(RniRw i>)(r) {R nl R n , v )(r') \(101' 0\L 0) | 2 (O) 

(A2) 

In order to c omp ute the elementary susceptibilities 
Xni ,n'V i n EqJA2), we have used the HF energies ex- 
tracted from ( |2.14 ) instead of the bare harmonic oscilla- 
tor ones. Further simplifications arise from the properties 
of the ClebschGordan coefficients and HO wave- functions 
R n i(r); in particular, explicit expressions for monopolar, 
dipolar and quadropolar excitations (L = 0, 1, 2) can be 
written. The monopolar free ph propagator reads 



^00 



( r y,n) = 



i 



(4tt) 5 



J2 (2l + l)R nl (r)R n , l (r) 



nln' 

Rai{r')Rn'i{r)xm,n'i{^) (A3) 

In this case we use the excitation operator O] = rf; this 
yields the following nonintcracting susceptibility 



A0 



= ~2 y"(2i + 1) [xm,n-nn(n + I + 1/2) 
at * — ' 

a nl 

i Xnl,n+ll (n+l)(n + l + 3/2)] (A4) 



which in the case of large enough number of particles, 
zero temperature and HO energies can be simplified to 
yield 

1 3 AT, 



(6JV ff ) 1/3 <S(fi- 2w a ) (A5) 



at 4 

Similarly, the dipolar propagator can be written as: 
G° o r\V) = j^tz E i(l + l)Rni(r)Rn'i+i(r) 

^ 77 ' nln' 

Rni(r')R (f2) + 1 R n i(r) 

Rnl-l{r)Rn I (r')R n i-i(r') Xnl,n> (A6) 

yielding a temperature independent free response 

1 3-ZVfy ( 1 1 



Xoi 



a^. 8n \ ft — uj a + ir/ ft + uj a + ir/ 



with a T = structure factor 
1 ZN a 



'-'Ol — 



(A7) 



(A8) 



[1] F. S. Cataliotti, E. A. Cornell, C. Fort, M. Inguscio, F. 

Marin, M. Predevelli, L. Ricci, and G. M. Tino, Phys. 

Rev. A 57, 1136 (1998); M. Prevedelli, F. S. Cataliotti, 

E. A. Cornell, J. R. Enscher, C. Fort, L. Ricci, G. M. 

Tino, and M. Inguscio, Phys. Rev. A 59, 886 (1999). 
[2] W. I. Mc Alexander, E. R. I. Abraham, N. W. M. Ritchie, 

C. I. Williams, H. T. C. Stoof, and R. G. Hulet, Phys. 

Rev. A 51, R871 (1995); M. O. Mewes, G. Ferrari, F. 

Schreck, A. Sinatra, and Ch. Salomon, Phys. Rev. A 61, 

R011403 (2000). 
[3] B. DeMarco, J. L. Bohn, J. P. Burke, Jr., M. Holland, 

and D. S. Jin, Phys. Rev. Lett. 82, 4208 (1999) 
[4] B. DeMarco and D. S. Jin, Science 285, 1703 (1999). 
[5] D. A. Butts and A. S. Roshsar, Phys. Rev. A 55, 4346 

(1997) . 

[6] J. Schneider and H. Wallis, Phys. Rev. A 57, 1253 (1998). 
[7] G. M. Bruun and K. Burnett, Phys. Rev. A 58,2427 

(1998) . 

[8] K. M0lmer, Phys. Rev. Lett. 80, 1804 (1998). 

[9] M. Amoruso, A. Minguzzi, S. Stringari, M. P. Tosi, and 

L. Vichi, Eur. Phys. J. D 4, 261 (1998). 
[10] L. Vichi, M. Inguscio, S. Stringari, and G. M. Tino, 

J. Phys. B: At. Mol. Opt. Phys. 31, L899-L907 (1998). 
[11] L. Vichi, M. Amoruso, A. Minguzzi, S. Stringari, and M. 

P. Tosi, Eur. Phys. J. D 11, 335-339 (2000). 
[12] L. Viverit, C. J. Pethick, and H. Smith, Phys. Rev. A 61 

053605 (2000). 
[13] A. Minguzzi and M. P. Tosi, Phys. Lett. A, in press. 



7 



[14] H. T. C. Stoof, M. Houbiers, C. A. Sackett, and R. G. 

Hulet, Phys. Rev. Lett. 76, 10 (1996). 
[15] M. Houbiers, R. Fewerda, H. T. C. Stoof, W. I. 

Mc Alexander, C. A. Sackett, and R. G. Hulet, Phys. Rev. 

A 56, 4864 (1997). 
[16] M. A. Baranov and D. S. Petrov, Phys. Rev. A 58, R801 

(1998) . 

[17] G. Bruun, Y. Castin, R. Dum, and K. Burnett, 

Eur. Phys. J. D 7, 433 (1999). 
[18] M. A. Baranov, JETP Lett. 70, 396 (1999). 
[19] W. Zhang, C. A. Sackett, and R. G. Hulet, Phys. Rev. A 

60, 504 (1999). 

[20] S. K. Yip and Tin-Lun Ho, Phys. Rev. A 59, 4653 (1999). 
[21] L. Vichi and S. Stringari, Phys. Rev. A 60, 4734 (1999). 
[22] G. M. Bruun and c. W. Clark, Phys. Rev. Lett. 83, 5415 

(1999) . 

[23] M. Amoruso, J. Meccoli, A. Minguzzi, and M. P. Tosi, 

Eur. Phys. J. D 7, 441 (1999). 
[24] M. Amoruso, J. Meccoli, A. Minguzzi, and M. P. Tosi, 

Eur. Phys. J. D 8, 361 (2000). 
[25] A. Fetter and D. Walecka, Quantum Theory of Many 

Particle Systems, McGraw Hill, New York, 1971. 
[26] G. Bertsch and S. F. Tsai, Phys. Rep. C 18, 126 (1975). 
[27] M. Casas and S. Stringari, J. Low Temp. Phys. 79, 135 

(1990). 

[28] M. M. Calbi and E. S. Hernandez, J. Low Temp. Phys. 

(2000) , in press. 

[29] E. S. Hernandez, J. Navarro, and A. Polls, Phys. Lett. B 
413, 1 (1997) and Nucl. Phys. A 627, 460 (1997). 

[30] S. M. Gatica, E. S. Hernandez, and J. Navarro, Phys. 
Rev. B 60, 15302 (1999). 

[31] D. Pines and P. Nozieres, The Theory of Quantum Liq- 
uids, Vol. II, Addison- Wesley, Readwood City, 1990. 

[32] W. Ketterle, D. S. Durfee, and D. M. Stamper-Kurn, in 
Bose-Einstein condensation of atomic gases, edited by 
M. Inguscio, S. Stringari, and C. E. Wieman, IOS press, 
Amsterdam, 1999. 

[33] W. Kohn, Phys. Rev. 123, 1242 (1961). 




FIG. 1. Monopolar dynamic structure factor (in arbitrary 
units) in the ETP case for JVi = 10 4 and uii = 2tt x 70s -1 . 
The upper plot is depicted in logarithmic scale and the lower 
one is a zoom of the main peak in linear scale. The thin line 
corresponds to the HF ph excitations, and thick and dashed 
lines indicate the RPA structure factor in the symmetric and 
antisymmetric channels, respectively. 




n/co 

FIG. 2. Dipolar dynamic structure factor (in logarithmic 
scale and arbitrary units) for the same conditions as in Fig.[j]. 
Thin, thick and dashed lines respectively correspond to the 
HF system, symmetric and antisymmetric channels. 
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FIG. 3. Monopolar density fluctuations for equal trap- 
ping potentials and populations (in arbitrary units) with 
Ni,2 = 10 4 . Solid and dashed lines respectively correspond to 
density and spin-like excitations at the poles of each channel. 
a^ 1 is the distance unit for the trapping potential of species 
1, a = (S/mui) « 0.53/im _1 . 
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FIG. 4. Dynamic structure factor for L = (in arbitrary 
units) for the symmetric channel of the interacting and free 
system with Ni = 10 4 and several concentrations, in thick 
and thin lines, respectively. Each column display a different 
range in the energy scale. 
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FIG. 5. Same as Fig.^ for the dipolar dynamic structure 
factor. 
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FIG. 6. Same as Fig.y for the quadrupolar dynamic struc- 
ture factor. 
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FIG. 7. Monopolar density fluctuations (in arbitrary units) 
for Ni = 10 4 . The upper, middle, and lower plots correspond 
to A = 3, 1,0.3 respectively, for <j\ = 7/2 (left column) and 
(72 = 9/2 (right column). Solid and dashed lines respectively 
indicate Sp^ and 5 pa ■ 
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FIG. 8. Same as Fig. |^ for dipolar excitations. 
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FIG. 9. Same as Fig. fj] for quadrupolar excitations. 

TABLE I. Shape parameters related to the HF probability 
density profiles P<r{r) (see text) for Ni = 10 4 and several 
concentrations. All distances are given in harmonic oscillator 
units a" 1 = (h/muji)^ 1 « 1.9^m. 
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CI 


ai 


cri 




cri 




-Pmax 


0.017 


0.022 


0.017 


0.018 


0.017 


0.015 


oS m „ 


5.64 


4.28 


5.64 


5.24 


5.64 


6.28 


C( -Radge 


8.71 


6.59 
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8.17 


8.79 


9.86 


FWHM 


5.04 


3.92 


5.04 


4.72 


4.96 


5.68 
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